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Abstract 

We study decoherence produced by a discrete environment on a charge Josephson qubit by 
introducing a model of an environment of bistable fluctuators. In particular we address the effect 
of 1// noise where memory effects play an important role. We perform a detailed investigation of 
various computation procedures (single shot measurements, repeated measurements) and discuss 
the problem of the information needed to characterize the effect of the environment. Although in 
general information beyond the power spectrum is needed, in many situations this results in the 
knowledge of only one more microscopic parameter of the environment. This allows to determine 
which degrees of freedom of the environment are effective sources of decoherence in each different 
physical situation considered. 

PACS numbers: 03.65.Yz, 03.67.Lx, 05.40.-a 
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A high degree of quantum coherence is crucial for operating quantum logic devices 
Solid state nanodevices seem particularly promising because of integrabilitv and flexibility 
in the design and several possible implementations have been proposed 0,3,0, 5,0]. Few 
recent experiments succeeded in detecting coherent dynamics in superconducting devices [3, 
y,y,[l^, but revealed limitations in the performances, due to decoherence. 

In a quantum logic device the interesting degrees of freedom are related to a given set 
of observables which we can prepare (write) and measure (read), and define the system. 
Their eigenstates | {<?«}) form the computational basis. The dynamics of a state \ip,t) = 
Xlgi qjsf ^gi--qNi'^) lily ■ ■ ^In) may be controlled if the Hamiltonian of the system is tunable. 
Loss of coherence is due to the fact that the Hilbert space of the device is much larger than 

nn 

the computational space [HI, 1121 ■ The additional degrees of freedom define the environment, 



which cannot be controlled and moreover little information is available on it. Decoherence 
is ultimately due to the entanglement between system and environment ^ , but in order to 
pursue this point of view one should be able study the full dynamics, in particular system- 
environment correlations. A less fundamental point of view, which we adopt here, is to 
associate decoherence to the loss of fidelity of a quantum gate. The environment blurs the 
output signal because it produces uncertainties in the phase relation between the amplitudes. 

The environment degrees of freedom may be "internal" to the device, which is a many- 
body object, or "external", belonging to the circuitry and to auxiliary devices. Internal 
decoherence is a serious problem in sohd state nanodevices, due to the presence of many low 
energy excitations. 

Investigation of the reduced qubit dynamics requires information about the environment. 
In several cases (weakly coupled environment [isj], harmonic oscillator environment 14, isl ]) 



the information is contained in the power spectrum of the environment operators coupled to 
the qubit 0- 

In this work we study a more general situation, an environment of quantum 
bistable fluctuators which may have memory on the time scale of the qubit dynamics, and 
may display the effect of non-gaussian correlations. A physical example are charged impuri- 
ties in substrates and oxides. These background charges (BC) produce for instance the 1/f 
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noise [17( observed in metallic tunnel junctions |18l I19I|. 

We will present here analytic results which elucidate several aspects of the BCs environ- 
ment. In particular we show that decoherence depends on details of the dynamics of the 
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environment beyond the power spectrum [20|. As a consequence it may differ for different 
gates, but in many cases the additional information required reduces to one microscopic 
parameter 21 1. 



MODEL FOR THE SYSTEM AND THE ENVIRONMENT 



Superconducting qubits mW are the only sol.d state .molenrentat.ons where coher- 
ence in a single qubit has been observed in the time domain [21 ISI iSlj and two-qubit systems 
are investigated 1^. Thus problems on decoherence can be posed in a realistic perspective. 
In all the above implementations it has been recognized that 1// noise is a major source 
of decoherence. We will focus on charge- Josephson qubits 0, where 1// noise produced by 
BCs trapped close to the device (Fig. We describe these BCs by introducing an impurity 
model. Many of our results and methods can be applied to different noise sources (eg. flux 
noise in flux- Josephson qubits 22, 13]) and other solid-state implementations (eg. gates 
based on the Coulomb interaction in semiconductor-spin qubits), but we will not discuss 
these topics here. 



A. The Superconducting Box 

The charge- Josephson qubit |5] is a superconducting island connected to a circuit via a 
Josephson junction and a capacitance C2 (Fig. (T]). The computational states are associated 
with the charge Q in the island. They are mixed by the Josephson tunneling. Level splittings 
can be tuned by the external voltage and under suitable conditions {Ec ^ Ej and 
low temperatures ksT <^ Ej ) only two charge states are important, which define the z 
component of a pseudospin. The the system implements a qubit with Hamiltonian 

'Hq = ^a,~^cT^ ; e{V,) = 4Ec{l - C^V.Ie) (1) 
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FIG. 1: A charge Josephson qubit in the presence of BCs located in the substrate or in the 
oxide close to the junction. Relevant scales are the charging energy Eq = e^/2(Ci + C2) and the 
Josephson energy Ej 

B. System Plus Environment Hamiltonian 



We describe each BC as a localized impurity level connected to a fermionic band 20Ll24, 



25|. For a single impurity the total Hamiltonian is 

^ = - I b'ba, + n' ; U' = e,b% + J^lTfeC^fe + h.c] + e.^c^ (2) 

k k 

Here Ti^ describes the BC alone: h (6^) destroys (creates) an electron in the localized level 
Ec and the electron may tunnel with amplitude to a band, described by the operators Cfc, 
c\ and the energies Sk- An important scale is the switching rate 7 = 2TiJ\f{t^\T\'^ [J\f is the 
density of states of the band, and jT^p ^ I^P)? which characterizes the relaxation regime 
of the BC. The BC determines a bistable extra bias v for the qubit, via the coupling term. 
For a set of BCs we generalize Eq. Q as follows 

i 

where extra bias operator is E{t) = Ylii'^iAbi- For simplicity the assumption that each 
localized level is connected to a distinct band has been made. 

C. Model for 1// Noise 

The environment in Eq. Q above is specified by the distribution of the switching rates 
7j which we choose in order to reproduce the 1// noise. The standard way is to assume 
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a distribution P(7) oc I/7 for 7 G [7m, 7a/]- Indeed if we look at the relaxation regime of the 
BC, the total extra polarization is a classical stochastic process E{t) with power spectrum 

Siuj) = J2s^{uj) ; S^ico) = IvUl-Sp) , (4) 

i 2 (7. +u}^) 

6p being the thermal average of the difference in the populations of the two states of the 
BC, the distribution ^(7) leads to to 1// noise, S{uj) = {7r(l —p'^)ndv'^/{4: In 10)} u^'^ for 
frequencies u G [7^, 7m] (^cZ is the number of fluctuators per noise decade). With the above 
choice of parameters, the Hamiltonian Eq. has been used to study decoherence due to 
1/f noise 



II. REDUCED DYNAMICS 



Our aim is to investigate the effect of the BC environment on the dynamics of the qubit. 
The standard road- map is to calculate the reduced density matrix [l4| of the qubit p^{t) = 
TTE{W{t)}, W{t) being the full density matrix. In the standard weak coupling approach 
a master equation for p'^(t) is written |13i], the environment entering its dynamics only via 
the power spectrum, Eq. (^. A master equation for p'^(t) can also be obtained by mode ling 
the environment by a set of harmonic oscillators. By using diagrammatic techniques [iJ, 
results of the standard weak coupling approach are obtained at lowest order in the 



couplings V, but it has been pointed out that higher orders are important for a 1// oscillator 



environment 



2G 
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The failure of the standard weak coupling approach is due to the fact that the 1/f envi- 
ronment includes fluctuators which are very slow on the time scale of the reduced dynamics, 
the approach being reliable only for BCs with Vi -C ji Q, • Rather than study higher 
orders in the perturbation series 2^, we prefer to use a different strategy. We study the 
models Eqs. fj2l3p by enlarging the system and considering only the bands as the environ- 
ment. The modified road-map consists in calculating in some approximation the reduced 
density matrix (RDM) p(t) = TTbands{W(t)} and then extract exactly p'^if). This allows to 
obtain results to all orders in the coupling v and to investigate details of various quantum 
gates. The disadvantage is that we have a larger system to deal with. We have investigated 
this problem with different techniques [2^, here we present a master equation approach, 
which covers many of the results we obtained before. 



A. Master Equation for a Single BC 



We split Eq. (j21) in system, TCq = Hq — | b^baz+Ecb'^b, and environment, He = ^kC^Ck, 
coupled by V = ^k[^kclb + h.c.]. The eigenstates (and eigenvalues) of Ho are product states 
of the form \qubit) \BC), namely \a) = \e+) |0) (-^^/2), \b) = {6-) |0) 
|c) = (-^]72 + ee), \d) = {n'/2 + ec). Here \e±) are the two 

eigenstates of an, the direction being specified by the polar angles 9 and = 0. Each set of 
qubit states corresponds to a value of 6^6 = 0, 1: the two level splittings are 



and Q' = ^/{e + vy + Ej, and finally cos^ = e/Q, sinO = Ej/Q, cosO' = {e + v)/Q', 
sin 9' = Ej/Vl'. 

In the basis of the eigenstates of Hq the master equation for the RDM in the Schrodinger 
representation reads 

dpij{t) 



dt 



mn Pmn it) (5) 



where uja is the difference of the eigenenergies and Rumn are the elements of the Redfield 



tensor 







-Sn, ^^^rnkir) e^^™- - 5.™ C< ^-.(t) c^ } (6) 
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The correlation functions are given by 

where iG^iu) = 7/(1 + e"'^'^), G'^iu) = G^^—u). Many of the coefficients Rijmn vanish. In 
particular the system of equations (0) splits in two blocks. The first contains the populations 
and the coherences pab and pd, together with their conjugate, i.e. the elements diagonal in 
the BC. The second block contains all the other coherences, which vanish identically if the 
initial RDM p(0) is diagonal in the BC, the physically relevant case for our purposes. 

In the standard secular approximation |l3| in the r.h.s. of Eq. © are retained only 
terms with coefficients Rjjjj, but this is not enough in our problem. To discuss possible 
approximations we first point out that we are interested in the case w -C fi, Q'. We have 
two important scales, namely Q' — Q and i7 ~ fi', the latter being much larger than the 
former. The secular approximation is valid if 7 ^ — fi, i.e. for an almost static BC. If 



'-f Q' — Q Q Q' we can no longer neglect the coefficients mixing pab and pcd (and 
those mixing the conjugates). We call this regime adiabatic. Finally in the fast BC regime, 
7 > Q', all the Rijmn must be present. Despite of this, the reduced dynamics of the qubit 



is just the result of lowest order approach 
anymore. 
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27|, so we will not discuss this regime 



B. Results for Adiabatic BCs 

In the adiabatic regime coherences are calculated using the following elements of the 
Redfield tensor 

Rabab = -^[i- - C^S - s'^S' + 1(0^10 + s'^w')] ] Rab cd = — + S - iw) 

Rcdcd = + c^^ + s^5' + i[c^w - s^w')] ; Rcdab = —{'^ - ^ - iw) 

where c = cos[(6' - 0')/2\, s = sm[{9 - 0')/2], S = t^a + tab, S' = t^a + tcb, w = Wca - Wdb, 
w' = Wda - Web and 

Here ipi^) is the digamma function. The coherences pab(t) and pcd(t) can be found in closed 
form by simply diagonalizing a 2 x 2 matrix. They allow to study the qubit coherence via 
{cTyit)) = —2Q[pabit) + Pcdit)]- The qubit coherence decays as exp{— r(t)}, where 

Pab{t) + Pcd{t) 



T{t) = - In 



(7) 



Pab{0) + Pcd{0) 

Here we present the analytic solution in a regime where the dynamics of the charge is 
not modified by the presence of the qubit, a case which is of interest for 1// noise. We find 

Pabit) + Pcdit) = e*(^+^^/2)* ^ {A(a)e-^(i-")* - e-^ (^+")*} (8) 

Za 



where a = y 1 — g'^ — 2ig6p — (1 — c^)(l — Sp^) and g = {Q' — ^) enter the decay rates 
whereas the prefactors are 

A{a) = {a + c^ - ig') pafe(O) + (a + + ig') pcd(O) ; g' = g + - c^) . 
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C. Two Regimes for the BC 



A rough analysis of Eq. (jS)) allows to draw a physical picture, li g ^ 1, a is substantially 
imaginary, reflecting the fact that a very slow BC mainly provides a static energy shift. 
Instead for (7 <^ 1 the decay rates acquire a real part thus a fast BC may determine an 
exponential reduction of the output signal. From a more quantitative analysis it emerges 
that BCs with g <^ 1, which we call weakly coupled, behave as a suitable set of harmonic 
oscillators with power spectrum given by Eq. (0)). They mainly produce the homogeneous 
broadening of the signal. Instead BCs with g ^ 1, which we call strongly coupled, give rise 
to memory effects and deviations of their statistical properties from those of an oscillator 
environment (cumulants higher than the second) are relevant [l^ . They mainly produce the 
inhomogeneous broadening of the signal. 

Finally we compare this result with other approaches. Notice flrst that Eq. (jHJ is valid 
if the power spectrum Si{u) (see Eq. (jU) is negligible at frequencies ~ f2. In this regime 
both the decay T{t) Eq. and the energy shift reproduce the results of the staiidard weak 
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coupling approach ^ if 7^ ^ 1. On the other hand, as in the work of Refs. 
may simulate the effect of the BC by a suitable set of quantum harmonic oscillators. The 
resulting decoherence rate for this oscillator environment is given by 
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(9) 



Jg=0 



where r(t) is given by Eq. (jHJ with thermal initial conditions for the BC. 



III. PURE DEPHASING 



For Ej = the environment only produces random fluctuations of the level splitting. 
This case, usually referred to as "pure dephasing", is special in that the Hamiltonian Q 
commutes with cr^. The charge in the island is then conserved and no relaxation occurs, but 
if we prepare the qubit in a superposition the system will dephase. 

If the initial density matrix is factorized, W{0) = we{0) ® p*^(0), it is possible to write 



an exact expression for the coherences only in terms of the environment 



pSW = P?i(0)e-^W-^^^W ; m = -ln|Tr^{^^(0)e^«-V^^^*}| 



(10) 
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where = ^-T^f - (r//2)E. This general expression may be further simphfied because 
individual charges contribute independently to T{t), T{t) = ^jrj(t), where Ti{t) is the 
dephasing due to a single charge resulting from the Hamiltonian Eq. Q. We notice that 
in this case 6 = 6', which implies that p%{t) = pab{t) + Pcd{t) and moreover the adiabatic 
approximation is exact, since the Redfield coefficients relating the above two coherences with 
all the other entries of the RDM vanish identically. Then Ti{t) is given by Eq. ((7j) and in 
the limit described by Eq. (jHl) we find the analytic form 

Tiit) = -\n A, e-* + (1 - A^) e"* (11) 



where A = [1 + (1 - igi)/a]poi + [1 + (1 + igi)/a]pii, gi = Vi/ji, at = ^1- gf - 2igi6pi, 
and poi (pii) is the probability that the i-th BC is initially empty (singly occupied). This 
result has been obtained using several different techniques in |2flL l2ll|. 

Finally by combining the equations of this section with Eq. Q we recover the exact 
result of a set of harmonic oscillators 

^ ,, f°°duj ^, , 1 — coscjt 

Toscit) = I — Siou) . (12) 







A. Single BC 

We study deviations of Eq. (fTT|) from the result for an oscillator environment, Eq. ()12j) . 
We consider different initial conditions for the BCs. Substantial deviations are clearly ob- 
served, except in the case of a weakly coupled BC {g = OA). In particular BCs with 
g > 1 induce slower dephasing compared to an oscillator environment with the same S{uj), 
a sort of saturation effect. Recurrences at times comparable with 1/v are visible in r(t). 
In addition strongly coupled charges show memory effects. We now study the effect of 
the initial conditions of the BC, expressed via 6pQ = Po — pi- In a single shot process 
5pQ = ±1 and r(t) describes dephasing during time evolution (homogeneous broaden- 
ing). Memory effects are apparent and in particular for 7^ <^ 1 the short time behavior 
is r(t) ^ vH'^{l - 6pl)/8 + jvH^ (1 + 26po6p - 35pg)/24 oc t^. A two level system is stiffer 
than a set of oscillators and indeed Tosdt) ~ v'^{l — 6p )/8t'^. On the other hand if we 
choose 6po = 6p in Eq. (fTT|). for very short times T(t) ~ Toscit)- This case corresponds 
physically to repeated measurements where the preparation of the BC is not controlled and 
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FIG. 2: Reduced T(t) due to a BC prepared in a stable state {6po = 1 solid lines and 5po = — 1 
dotted lines) and in a thermal mixture {Spo = 6p dashed lines) for the indicated values of g = v/^. 
Inset: longer time behavior for stable state preparation. The curves are normalized in such a way 
that the oscillator approximation for all of them coincides (thick dashed line) 

determines slightly different characteristic frequency for the qubit. This sort of inhomoge- 
neous broadening adds to decoherence during the individual time evolutions and determines 
a faster decay of Tit). 

To summarize a weakly coupled BC, (? ^ 1, behaves as a source of gaussian noise, 
decoherence depending only on the power spectrum of the fluctuator, whereas decoherence 
due a strongly coupled BC, g ^ 1, displays saturation effects and dependence on the initial 
conditions of the BC. 

B. 1// Noise in Single Shot Measurements 

The set of BCs producing 1/f noise contains both weakly and strongly coupled fiuctua- 
tors, so no typical time scale is present. A very large number of slow fiuctuators is present, 
and it is not a priori clear how saturation manifests. In Fig. El we show the results for a 
realistic sample. We choose initial conditions 6pQj = ±1 randomly distributed on the set 
of BCs with (l/N) ^Poj ~ ^P- We checked that different microscopic realizations of 
this conditions give roughly the same total initial extra bias E{0) and the same dephasing. 
We calculate dephasing during the time evolution, i.e. the average signal of several single- 
shot experiments, where E[0) is re-calibrated before each experiment. This ideal protocol 
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FIG. 3: r(t) due to 1// noise in the range uj E [7^, 10^^ Hz], for decreasing 7^ (solid lines, the label 
is the number of decades included). Slower BCs saturate whereas this not happens for an oscillator 
environment (dashed lines). Parameters {y = 9.2 x IQ^ Hz, rid = 1000) give the experimental noise 
levels and reproduce the observed decay of the echo signal ly] . Couplings Vi are distributed with 
Af /R = 0.2 



minimizes the effects of the environment. 

We now perform a spectral analysis of the effects of the environment by adding slower 
fiuctuators decade after decade, in a way such that 1// noise with the same amplitude A is 
present for uj G [7m, 10^^ Hz\. We see in this example that dephasing is given by BCs with 
7j > lO^Hz ^ I'^I/IO. The overall effect of the strongly coupled BCs (7^- < Itil/lO) is minimal, 
despite of their large number, thus low-frequency noise saturates. Instead Toscit), Eq. (fT^ 
does not saturate at low frequencies. We notice finally that for this protocol T(t) is roughly 
given by Toscit) provided low frequencies are cut off at u ~ l^;!. Thus dephasing depends 
essentially on a single additional parameter besides the power spectrum S{uj), namely the 
average coupling |t>| or equivalently the number of charges producing a decade of noise, n^. 

Our results are not very sensitive to the value of 11^ we choose. Indeed for constant 

amplitude A we must keep constant "^^vf ~ n^f^, so the effective low-frequency cutoff 
1 — r —1/2 

~ |f I varies as . For ^ 00 the low-frequency cutoff goes to zero, as in the result Eq. 
(fT^ for the oscillator environment. 
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C. Repeated Measurements and Inhomogeneous Broadening 



Single-shot measurements are a goal for experimental research, but presently available 
protocols involve repeated measurements. For instance in [3, the time evolution proce- 
dure is repeated ~ 10^ times and the total current due to the possible presence of the extra 
Cooper pair in all the repetitions is measured. The signal is then the sum over different 
possible time evolutions of the BCs with initial conditions which are also randomly fluctu- 
ating. This additional blurring of the output signal results in a faster decay for T{t). As 



explained in 



2l|, this is accounted for by letting in r(t) (i.e. in Eqs. (fTT|) ) 6poi = Sp^^tr 
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FIG. 4: Different averages over Spoj for 1// spectrum reproduce the effect of repeated measure- 
ments. They are obtained by neglecting (dotted lines) or account for (solid lines) the strongly 
correlated dynamics of 1// noise. The noise level of |l9| is used, by setting \v\ = 9.2 x W^Hz, 
= 10^, 7m = IHz, 7m = W^Hz. Dashed lines are the oscillator approximation with a lower 
cutoff at uj = min{|w|, 1/tm} 



the average of the values of 6pj sampled at regular times ta for the overall measurement 
time tm- As a rough estimate we let Sp^itm) = ±1 if Ti^m < 1 and Sp^itm) = Sp^ for 
7jtm > 1, and consider only the case of long overall measurement time |f |tm ^ 1- In this 
case BCs with 7 < are saturated and are not effective, whereas for the other BCs, 

being averaged, we may take F(^)(t) ^ FE(t) for small enough times. The result would be 
F(t) ~ Jy^ doj S{uj){l — cos cijt) /(ttcj^) and would proof the recipe proposed by Cottet et 

al. y . 
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In Fig. m (dotted lines) we show that indeed dephasing calculated as outlined above, 
is roughly given at short times by the oscillator environment approximation with a lower 
cutoff taken at uj ^ min{|f |, 1/tm}- We also show results with a different averaging procedure 
Sp^{tm) = ^/tm /o™ dt6p{t) which takes into account the strongly correlated dynamics of 1// 
noise (solid lines in Fig. IH). These correlations do not affect the results except possibly for 



D. Charge Echo 

Echo-type techniques have been recently suggested 0,13 and experimentally tested [3| 
as a tool to reduce inhomogeneous broadening due to the low-frequency fiuctuators of the 
1/f spectrum. In the experiment of 19] the echo protocol consists of a 7r/2 preparation 
pulse, a TT swap pulse and a 7r/2 measurement pulse. Each pulse is separated by the delay 
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time t. We calculated the decay of the echo signal using a semi-classical approach 
This result allowed to estimate the parameters we have used in this work by comparing with 

Q. 

The decay depends very weakly on initial conditions of the BCs and is well reproduced by 
the oscillator environment approximation. This means that in the experiment [l^ the echo 
procedure actually cancels the effect of strongly coupled charges, this conclusion being valid 
as long as the delay time is short, t\v\ <^ 1. We remark that in the regime of parameters 
we consider, for given noise amplitude the echo signal is strongly dependent on the high 
frequency cutoff and a detailed analysis may give information on the actual existence of 
BCs switching at rates comparable with ~ 10 GHz. 



IV. DECOHERENCE FOR A GENERIC WORKING POINT 



The most effective strategy for defeating 1/f noise has been implemented in the exper- 



iment by Vion et al. 



. It is useful to explain it in a non pictorial way by considering 



decoherence as given by the weak coupling approach 
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ro(t) 



1 



s{o) cos' e + - s{n) sin' e 



t. 



(13) 



.2 ^ ' 4 

Even if this formula does not hold for 1/f noise, it indicates that the dangerous adiabatic 
term containing 5(0) may be eliminated (in lowest order) if one operates at 6* = 7r/2. 
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Small deviations from this point may determine a dramatic increase of decoherence j8|. 
Understanding decoherence in a generic operating point may help in designing more flexible 
gates, or to implement different strategies as computation with geometric phases 0]. 



A. Single BC 

First we consider Eq. (jH)) and we notice that by changing the working point 6', a strongly 
coupled BC may become weakly coupled and vice-versa. For simplicity we let KT ^ Ec-, 
then a ~ a/c^ — g"^. The BC is weakly coupled if 7 7c = {Vt' — Vl)/c^. Thus, the threshold 
value decreases if we go from 6 = Q to 6 = -n /2 (see Fig. IIV|1 . 

To get some insight in the problem of decoherence away from 6 = 71 /2 we consider r{t) 
in the adiabatic regime Eqs. ()7I8|) . Results plotted in Fig. IHl show r(t) parametrized by 
6 for four values of the parameter f/7. For reference we plot also (dashed lines) ro(t), 
Eq. The top left panel shows a weakly coupled BC for which r{t) roughly follows 

ro(t). The other panels show BCs with f > 7 which turn from strongly coupled to weakly 
coupled increasing 9 G [0, 7r/2]. For these latter BCs saturation is less effective in suppressing 
decoherence when the operation point is close to 6 = 7r/2, which may be an indication of 
the fact that their effect is more sensitive to deviations from the optimal point. 




0.2 0.4 0.6 0.8 1 1.2 1.4 



e 

FIG. 5: The threshold value for a BC behaving as weakly coupled depends on the operating point 
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FIG. 6: Reduced T{t) (solid lines) plotted for various f/7 = 0.05, 1.0, 2.0, 4.0. In each panel 
curves correspond, from top to bottom, to 9 = tt/2, 0.45 vr, vr/S, 7r/4, 0. For comparison TQ{t) is also 
plotted (dashed lines). In the units used, ro(t) is the same in each panel 

B. 1/f Noise at the Optimal Point 

For a set of BCs the road-map for the reduced dynamics outlined in Sec. IHlis not easily 
implemented by generalizing the master equation Eq. (0). In we used the Heisenberg 
equations of motion, where we factorize correlations between the qubit and the bands, and 
between different BCs. We are left with a set of 3(A^ + 1) coupled differential equations, 
where N is the total number of BCs. These approximations are expected to work for small 
Vi but in effect they give accurate results for general values of Qi even if Vi/Ej is not very 
small, as we checked by comparing with numerical evaluation of the master equation for the 
qubit with one and two BCs. We study decoherence at the optimal point 6 = 7i/2 {Q = Ej) 
via the Fourier transform of {crz{t)). We consider a set of weakly coupled BCs in the range 
[10~^,10]-Ej which determines 1/f noise around the operating frequency with amplitude 
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(b) 




FIG. 7: (a) The Fourier transform crz{uj) for a set of weakly coupled BCs plus a single strongly 
coupled BC (solid line). The separate effect of the coupled slow EC alone {go = 8.3, dashed line) 
and of the set of weakly coupled BCs (dotted line), is shown for comparison. In the inset the 
corresponding power spectra: notice that at a; = Ej the power spectrum of the extra charge alone 
(dashed line) is very small. In all cases the noise level at Ej is fixed to the value S{Ej)/Ej ~ 
3.18 X 10^"^. (b) The Fourier transform az{u)) for a set of weakly coupled BCs plus a strongly 
coupled BC {vq/'Jq = 61.25) prepared in the ground (dotted line) or in the excited state (thick 
line) 

of the typical measured spectra Q, [l^ extrapolated at GHz frequencies. Dephasing due 
to this set of BCs agrees with the weak coupling result Eq. ()13|). Then we add a slower 
(and strongly coupled go = Vq/jq = 8.3) BC, which should produce no effect according to 
Eq. ((131), Fig- Ht- We find that the strongly coupled BC alone determines a dephasing rate 
comparable with that of the weakly coupled BCs and the overall dephasing rate is more than 
doubled. This result shows that slower charges 7 ^ play a role in dephasing. Moreover 
information beyond S{uj) is needed, as we checked by showing that sets of charges with 
different and Vi but the same S{uj) yield substantially different values of the decoherence 
rate. Decoherence is larger if BCs with g>l are present in the set. In summary these 
results show that Eq. (fT!?|) underestimates the effect of strongly coupled BCs and also that 
decoherence at the optimal point may be substantial even if the 1// spectrum does not 
extend up to frequencies ~ Ej. 

If we further slow down the added BC we find that increases toward values ~ 70, 
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the switching rate of the BC. This indicates that the effect of strongly coupled BCs on 
decoherence tends to saturate, in analogy with the results for pure dephasing. In this 
regime, we observe also memory effects related to the initial preparation of the strongly 
coupled BC (see Fig. Ub)- Again we expect a dependence of dephasing on the protocol of 
the quantum gate. 

V. CONCLUSIONS 

In conclusion we have studied dephasing due to charge fluctuations in solid state qubits. 
For a fluctuator environment with 1// spectrum memory effects and higher order moments 
are important so additional information on the environment is needed to estimate dephasing. 
In this case the additional information of the environment needed depends on the protocol 
but often reduces to a single parameter. A new energy scale emerges, the average coupling 
|f I of the qubit with the BCs, which is the additional information needed to discuss single 
shot experiments (alternatively one should know the order of magnitude of Ud). For repeated 
experiments the relevant scale is instead min{|f |, 1/tm} where is the overall measurement 
time. We point out that these scales emerge directly from the study of the dynamics of the 
model, and not from further assumptions. Finally echo measurements are sensitive to the 
high-frequency cutoff •jm of the 1// spectrum. 

It is interesting to notice that some result relative to an environment of quantum harmonic 
oscillators for arbitrary qubit-environment coupling may be obtained as a limit of the discrete 
environment in the semi-classical regime, and depend on the classical statistical properties of 
an equivalent random process with no reference to the quantum nature of the environment. 

Our results are directly applicable to other implementations of solid state qubits. Joseph- 
son flux qubits Isoj] suffer from similar 1// noise, originated from trapped vortices. Our 
model applies if cx^ represents the flux. Also the the parametric effect of 1// noise on the 
coupling energy of a Josephson junction 23] can be analyzed within our model, as long 
as individual fiuctuators do not determine large variations of Ej. In this case it may be 
possible that the same sources generate both charge noise and fluctuations of Ej. This can 
be accounted for by choosing the "noise axis" as the z axis. 

An important issue is to understand dephasing near the optimal operating points of the 
qubit P . Low- frequency noise can also be minimized by echo techniques, but the flexibility 
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in the implementation of gates is greatly reduced. A possibility is to implement quantum 
computation using Berry phases p, where the design of gates includes an echo procedure. 
Frequency shifts can be calculated within our model but a reliable analysis of the effect of 
1// noise is still missing. 

Finally we mention that the sensitivity of coherent devices may be used to investigate high 
frequency noise 0|. In particular an accurate matching between measured inhomogeneous 
broadening, echo signal and relaxation may give reliable information on the actual existence 
of BCs at GHz and on the high-frequency cutoff of the 1// spectrum. 
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